Effect of intrinsic spin relaxation on the spin-dependent cotunneling transport 

through quantum dots 
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Spin-polarized transport through quantum dots is analyzed theoretically in the cotunneling 
regime. It is shown that the zero-bias anomaly, found recently in the antiparallel configuration, 
can also exist in the case when one electrode is magnetic while the other one is nonmagnetic. Phys- 
ical mechanism of the anomaly is also discussed. It is demonstrated that intrinsic spin relaxation in 
the dot has a significant influence on the zero-bias maximum in the differential conductance - the 
anomaly becomes enhanced by weak spin-flip scattering in the dot and then disappears in the limit 
of fast spin relaxation. Apart from this, inverse tunnel magnetoresistance has been found in the 
limit of fast intrinsic spin relaxation in the dot. The diode-like behavior of transport characteristics 
in the cotunneling regime has been found in the case of quantum dots asymmetrically coupled to 
the leads. This behavior may be enhanced by spin-flip relaxation processes. 



PACS numbers: 72.25.Mk, 73.63.Kv, 85.75.-d, 73.23.Hk 



I. INTRODUCTION 



Electronic transport through nano-scale systems cou- 
pled to ferromagnetic leads reveals new phenomena which 
arise due to the interplay of charge and spin degrees of 
freedom. In the case of metallic grains or large quantum 
dots, this leads to magnetic-dependent Coulomb block- 
ade and Coulomb oscillations phenomena 1 *^, spin accu- 
mulation, negative differential conductance, and inverse 
tunnel magnetoresistance (TMR) effecli 4 *^. Transport 
through a quantum dot with a few electrons reveals fur- 
ther features resulting from discreteness of the dot energy 
spectrum and from electron correlationsZ&SiiSiii. When 
the dot is weakly coupled to the leads, one can distinguish 
two different regimes of electronic transport: sequen- 
tial and cotunneling ones 12,13,14 . Sequential transport 
through a single-level quantum dot coupled to ferromag- 
netic leads was studied quite recently for collinea r 15 : 16 
and non-collineaniLiS* 1 ^ configurations of the magnetic 
moments of electrodes. Spin polarized transport in the 
cotunneling regime has also been addressed22i2ii2£. It has 
been shownSl that the interplay of various spin depen- 
dent cotunneling processes gives rise to anomalous behav- 
ior of the differential conductance when the leads' mag- 
netic moments are antiparallel. This, in turn, leads to 
the corresponding zero-bias anomaly in the tunnel mag- 
netoresistance. Moreover, the even-odd electron number 
parity effect in tunnel magnetoresistance has also been 
founds. In both cases intrinsic spin relaxation in the 
dot was neglected, and spin in the dot could relax only 
due to cotunneling events. 

In this paper we consider transport through single-level 
quantum dots coupled to ferromagnetic leads in the co- 
tunneling regime. In particular, we show how the in- 
trinsic spin relaxation in the dot affects the dot occupa- 
tion probabilities, spin accumulation, conductance and 
the TMR effect. We analyze the cases of quantum dots 



with degenerate and non-degenerate energy levels, which 
are symmetrically (equal spin polarizations of the leads) 
and asymmetrically (unequal spin polarizations of the 
leads) coupled to the leads. The latter case is of par- 
ticular interest as the corresponding transport charac- 
teristics are highly asymmetric with respect to the bias 
reversal. We show that this diode-like behavior may 
be enhanced by intrinsic spin-flip relaxation processes in 
the dot. Moreover, we show that the zero-bias anomaly, 
found previously^ in the antiparallel configuration, can 
occur in any magnetically asymmetric system, in partic- 
ular when one electrode is nonmagnetic while the other 
one is ferromagnetic. 

The systems considered in this paper may be re- 
alized experimentally in various ways, including self- 
assembled dots in ferromagnetic semiconductors^, ultra- 
small metallic particles^, granular structures^, carbon 
nanotubes^^&Sii, and other molecules^ .. In the case of 
quantum dots coupled to nonmagnetic leads, differential 
conductance in the cotunneling regime has also been mea- 
sured experimentally 31 . 

In section 2 we present the model and theoretical for- 
mulation of the problem. Numerical results for symmet- 
ric and asymmetric systems are shown and discussed in 
sections 3 and 4, respectively. Final conclusions are given 
in section 5. 



II. MODEL AND METHOD 

We consider a single-level quantum dot coupled to 
ferromagnetic leads and for simplicity restrict our con- 
siderations to the case when magnetic moments of the 
leads are either parallel or antiparallel. The system is 
described by the Hamiltonian H = + + Hd + 
Ht. The first two terms represent the left and right 
reservoirs in the noninteracting particle approximation, 
H v = Y,kcr £ vkacl kcr Cv k(7 , where e vka denotes energy 
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of an electron with the wave number k and spin a 
in the lead v — L, R. The dot Hamiltonian H-q in- 
volves two terms, = J2a=i \ e ad\d a + Ud^d^d^di, 
where the first term describes noninteracting electrons 
in the dot, whereas the second one takes into account 
Coulomb interaction with U denoting the correlation 
parameter. The tunneling Hamiltonian reads, ifr = 

E^L.itEfea {Tycl ka d a +T„*4<w), where T v denotes 
the corresponding tunnel matrix elements. Coupling of 
the dot to external leads can be described by the parame- 
ters VI = 2it\Tv\ 2 pvc, with p va being the spin-dependent 
density of states in the lead v. When defining the spin po- 
larization of lead v as p„ = (r+— r~)/(r++r~), the cou- 
pling parameters can be written as Tt ' = r„(l ±p„), 
with = (r+ + r-)/2. Here, T+ and T~ correspond to 
spin-majority and spin-minority electrons, respectively. 
In the following we assume IY = Tr =T/2. In the weak 
coupling regime, typical values of the dot-lead coupling 
strength V are of the order of tens of peV, whereas the 
measurements are usually carried out at temperatures of 
the order of tens of mK 31 . 

In order to analyze electronic transport in the cotun- 
neling regime we employ the second-order perturbation 
theory 1 ?' 3 -. The cotunneling rate for electron transition 
from the left to right leads is then given by 
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where &j and e/ denote the energies of initial and fi- 
nal states, is the state with an electron in the 
lead v, whereas is a virtual state with s v denot- 
ing the corresponding energy. Generally, one can distin- 
guish between single-barrier and double-barrier cotunnel- 
ing as well as between cotunneling processes that leave 
unchanged magnetic state of the dot and those which re- 
verse spin of the dot. The non-spin- flip cotunneling pro- 
cesses are elastic and coherent, whereas the spin- flip ones 
are incoherent (and inelastic when the dot level is spin 
split). Unlike the double-barrier processes, the single- 
barrier ones do not contribute directly to electric current. 
However, they can change the occupation probabilities, 
and this way also the current flowing through the system. 

In numerical calculations we have included all possi- 
ble cotunneling events. Furthermore, we have allowed 
for intrinsic spin relaxation processes in the doti^. Such 
processes can result, for instance, from spin-orbit inter- 
action in the dot or coupling of the electron spin to nu- 
clear spins. We will not consider a particular microscopic 
mechanism of the intrinsic spin- flip processes, but simply 
assume their presence and describe them by the relevant 
spin- flip relaxation time r s f. The relaxation processes 
have been then taken into account via a relaxation term 
in the appropriate master equation for the occupation 
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where (3 = 1/(/cbT 1 ), P a denotes the probability that the 
dot is occupied by a spin-er electron, and 7^ t 7 tT is the 
cotunneling rate from lead v to lead v' with a change 
of the dot spin from a to a (a = —&). The last term 
describes the spin relaxation processes. In the case of 
spin-degenerate dot level this term is reduced to — (P a — 

We consider the case when the dot is singly occupied 
at equilibrium (e a < 0, e a + U > 0) and the system is in 
a deep Coulomb blockade regime; T, k^T <C |e CT |,£ CT + U. 
In such a case the sequential tunneling is exponentially 
suppressed, and cotunneling gives the dominant contri- 
bution to electric current. Both non-spin-flip and spin- 
flip cotunneling processes are then allowed and have to 
be taken into account. 

In the following we will distinguish between the fast 
and slow spin relaxation limits. The former (latter) limit 
corresponds to the situation when the time between suc- 
cessive cotunneling events, r cot , is significantly longer 
(shorter) than the intrinsic spin relaxation time r s f. A 
typical spin relaxation time for quantum dots can be 
relatively long, up to fisM^l. On the other hand, the 
time between successive cotunneling events can be esti- 
mated taking into account the fact that the rate of spin- 
flip cotunneling is generally larger than that of non-spin- 
flip cotunneling (for a finite parameter U). Assuming 
e-t — e — s. one then finds 
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with A — max{|ey|, knT} and h — 2nh. Assuming typi- 
cal parameters^, one can roughly estimate r cot to range 
from 10 -3 ns to 1 ns. Thus, having spin relaxation time 
in the dot and the system parameters entering Eq. J2J, 
one can estimate whether the system's behavior corre- 
sponds to the fast or slow spin relaxation regimes. Both 
the time between successive cotunneling events and the 
spin relaxation time depend on the system parameters 
and can be tuned, for example by applying a weak mag- 
netic field, varying coupling to nuclear spins, etc. Thus, 
the regimes of fast and slow spin relaxation are achievable 
experimentally. 

Electric current flowing through the system from the 
left to right leads is given by 
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In the following we discuss the influence of intrinsic spin 
relaxation on differential conductance and tunnel mag- 
netoresistance of the quantum dot in the cotunneling 
regime both in the presence and absence of external mag- 
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netic field B. The TMR effect is defined as 

_ ^AP 



TMR = 



IaP 



(5) 



with Ip (Zap) denoting the current flowing through the 
system in the parallel (antiparallel) magnetic configura- 
tion at a constant bias voltage V applied to the system. 



III. QUANTUM DOTS SYMMETRICALLY 
COUPLED TO THE LEADS 



First, let us proceed with the case of quantum dots 
coupled to ferromagnetic leads with equal spin polariza- 
tions = pr = p) . Our considerations will also include 
the cases of quantum dots with non-degenerate and de- 
generate energy levels. The former situation is discussed 
in the following subsection. 



A. The case of spin-degenerate dot level 

In the case of B = the dot level is spin degenerate, 
e-f = £j_ = e. The differential conductance G = dl/dV 
for parallel and antiparallel magnetic configurations is 
shown in Fig. da) for the case of no intrinsic spin relax- 
ation in the dot (l/r s f — > 0). In the parallel configuration 
one finds typical parabolic behavior of the conductance 
with increasing transport voltage. However, this is not 
the case for antiparallel configuration, where a maximum 
in the differential conductance appears in the small bias 
voltage regime. This anomalous behavior of the differen- 
tial conductance, in turn, leads to respective minimum 
in the tunnel magnetoresistance, as shown by the solid 
curve in Fig. Pfo). 

Physical mechanism of the zero-bias anomaly was con- 
sidered in Ref. |2(| for the case of no intrinsic spin re- 
laxation. It was pointed out that the anomaly results 
from the relative difference between the fastest cotunnel- 
ing events (which involve only the majority electrons of 
the leads) contributing directly to the charge current in 
the parallel and antiparallel configurations. In the par- 
allel configuration, the main contribution to electric cur- 
rent comes from the non-spin-flip cotunneling processes, 
whereas spin-flip cotunneling is dominant for antiparallel 
alignment. In addition, the spin asymmetry in tunneling 
processes leads to a nonequilibrium spin accumulation in 
the dot in the antiparallel configuration, which increases 
with increasing bias voltage. The spin accumulation, de- 
fined as (P-f — Pi)/ 2, is shown in Fig.|21as a function of the 
bias voltage for different spin relaxation times. For the 
parameters assumed here and for negative bias (eV < 0, 
current flows from right to left, while electrons flow from 
left to right), there is a larger probability to find in the 
dot a spin-up electron than a spin-down one. This, in 
turn, diminishes the contribution coming from the fastest 
spin-flip cotunneling processes. However, there are also 
single-barrier spin-flip cotunneling events which reverse 
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FIG. 1: The differential conductance in the parallel and an- 
tiparallel configurations (a,b) and tunnel magnetoresistance 
(c) as a function of the bias voltage for different spin re- 
laxation r — h/(T s fT). The parameters are: ksT — 0.2T, 
e = -15r, U = 30r, and p L = m = 0.5. 



spin of the dot and open the system for the fastest trans- 
port processes. The rate of single-barrier processes is 
proportional to temperature (but independent of the ap- 
plied voltage), whereas that of double-barrier processes 
increases with the transport voltage, being roughly inde- 
pendent of temperature. Thus, when \eV\ < k^T, the 
single-barrier processes play a significant role - they de- 
crease the spin accumulation in the dot and open the 
system for the fastest cotunneling events. Consequently, 
the main contribution to current comes then from the 
most probable spin- flip cotunneling processes. When 
\eV\ 3> k-gT, the relative role of single-barrier processes 
is diminished due to an increased role of double-barrier 
cotunneling, and the differential conductance decreases. 
The zero-bias maximum in differential conductance ap- 
pears then as a result of the competition between single- 
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be then expressed as 
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FIG. 2: The spin accumulation in the antiparallel config- 
uration as a function of the bias voltage for different spin 
relaxation r = /i/(r s fT). The parameters are the same as in 
Fig. EI 



and double-barrier cotunneling processes. 

The anomalous behavior of differential conductance 
in the antiparallel configuration results from spin asym- 
metry of tunneling processes. This asymmetry leads 
to nonequivalent occupation of the dot by spin-up and 
spin-down electrons (spin accumulation), as illustrated in 
Fig. From the above discussion follows that the spin 
accumulation is crucial for the occurrence of the zero- 
bias anomaly. Thus, one may expect that intrinsic spin- 
flip processes in the dot should suppress the anomaly. 
Indeed, intrinsic spin relaxation in the dot significantly 
modifies the conductance maximum at zero bias, and 
totally suppresses this anomalous behavior in the limit 
of fast spin relaxation. This behavior is displayed in 
Fig. ^b) , where different curves correspond to different 
values of the parameter r defined as r = h/{r s fT). Thus, 
r = describes the case with no intrinsic spin relaxation, 
whereas the curves corresponding to nonzero r describe 
the influence of intrinsic relaxation processes. First of 
all, one can note that small amount of intrinsic spin-flip 
processes enhances the zero-bias anomaly [see the curve 
for r — 10~ 2 in Fig. ^b)]. This is because such pro- 
cesses play then a role similar to that of single-barrier 
spin-flip cotunneling. In the case of fast spin relaxation, 
on the other hand, the spin accumulation is suppressed 
and the anomaly disappears, as can be seen in Fig. EJb) 
and Fig. [21 for r = 1 . It is also worth noting that spin- 
flip processes in the dot enhance the overall conductance 
in the antiparallel configuration. In the parallel config- 
uration, however, the differential conductance does not 
depend on intrinsic relaxation. This is because there is 
no spin accumulation in the parallel configuration, and 
consequently the intrinsic spin-flip processes in the dot 
do not play any role. 

In the case of a deep Coulomb blockade one can take 
only the lowest order corrections in x/y, with x — 
ksT, \eV\ and y = \e\,e + U, and derive some approx- 
imate formulas for the differential conductance. For the 
parallel configuration, the differential conductance can 
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In the antiparallel configuration and for \eV\ <C k B T, the 
conductance is given by the expression 
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which describes the maximum of the differential conduc- 
tance at zero bias for arbitrary value of r. On the other 
hand, when \eV\ 3> k B T and r<l one finds the formula 
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which in turn approximates the local minimum value of 
the differential conductance. The expression for G m ?, 
for arbitrary spin relaxation time is too complex to be 
presented here. It is clearly evident from the above ex- 
pressions that the conductance depends on intrinsic spin 
relaxation only in the antiparallel configuration, in agree- 
ment with numerical results. The variation of the linear 
conductance in the antiparallel configuration with the pa- 
rameter r is displayed in Fig. Eta) for several values of 
temperature. As one can see, there is a clear crossover 
between the limits of fast and slow spin relaxation in the 
dot. The crossover depends on temperature and is shifted 
towards shorter relaxation times as the temperature in- 
creases. 

In order to describe the zero-bias anomaly in the dif- 
ferential conductance in the antiparallel configuration, it 
is useful to introduce the relative height of the maximum 
at V = 0, defined as x G = (G^ax ~ G m f n )/G m f n . The 
variation of xq with the spin relaxation time is shown in 
Fig. Etb) for different temperatures. First of all, one can 
see that the relative height depends on r in a nonmono- 
tonic way. In the limit of slow spin relaxation and low 
temperature one finds xq — Ap 2 / (3 — p 2 ), whereas in the 
limit of fast spin relaxation xq tends to zero, indepen- 
dently of temperature. Indeed, in the limit of fast spin 
relaxation G m a X an( i G^?, are equal and given by 
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On the other hand, the relative height exhibits a maxi- 
mum for a certain value of spin relaxation. This maxi- 
mum, in turn, is decreased and shifts towards shorter re- 
laxation times with increasing temperature, see Fig.^b). 
One can show that the maximum occurs when the relax- 
ation time is approximately given by 

(10) 
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When comparing Eqs 10 and l|10|) one can conclude 
that the maximum in the relative height of the conduc- 
tance anomaly occurs when the spin relaxation time r 8 f is 
roughly equal to the time between successive cotunneling 
events, r cot , in the zero bias limit. This also means that 
the maximum separates the fast and slow spin relaxation 
regimes, in agreement with Fig. |2| The enhancement of 
xg can range from a few percent for small spin polariza- 
tion of the leads, p < 0.2, to several hundred percent for 
p > 0.9. For the parameters assumed in Fig. Bib), the 
relative height of the zero-bias maximum in differential 
conductance for t s { = T^ ax is enhanced by about 65% as 
compared to the case of r s f — * oo. Assuming typical pa- 
rameters for the quantum do1*2I, one can estimate r s ™ ax to 
be T s ™ ax sw 50ns, which is accessible experimentallyi2&2i 

Intrinsic spin-flip scattering in the dot modifies con- 
ductance in the antiparallel configuration, and thus also 
the tunnel magnetoresistance, as shown in Fig. ^c). 
Since the zero-bias maximum in conductance is sup- 
pressed in the fast spin relaxation limit, the correspond- 
ing dip in TMR at small voltages is also suppressed by 
the intrinsic relaxation processes. More specifically, the 
dip in TMR broadens with increasing r and disappears 
in the limit of fast relaxation (see the curve for r = 1). 

An interesting and new feature of TMR in the pres- 
ence of spin-flip scattering in the dot is the crossover 
from positive to negative values when r increases, as il- 
lustrated in Fig. ^c). Thus, the difference between con- 
ductances in the parallel and antiparallel magnetic con- 
figurations persists even for fast spin relaxation in the 
dot, contrary to the sequential tunneling regime, where 
such a difference disappears^. This seemingly counter- 
intuitive behavior can be understood by taking into ac- 
count the following two facts: (i) absence of spin accu- 
mulation in the dot for fast spin relaxation (Pf = Pj), 
and (ii) difference in the fastest cotunneling processes 
contributing to the current in the two magnetic config- 
urations. As the point (i) is rather obvious in view of 
the results presented in Fig. |2 the point (ii) requires 
additional clarification. The fastest double-barrier co- 
tunneling processes involve only the majority-spin elec- 
trons of the two leads - thus, in the parallel configura- 
tion the fastest cotunneling processes are the non-spin- 
flip ones. They take place either via the empty-dot vir- 
tual state (for one orientation of the dot spin) or via the 
doubly occupied dot virtual state (for the second orien- 
tation of the dot spin). The dominant contribution to 
the current is then proportional to 1/e 2 + l/(e + U) 2 . 
On the other hand, in the antiparallel magnetic configu- 
ration the fastest cotunneling processes are the spin-flip 
ones, which can occur only for one particular orienta- 
tion of the dot spin. However, for this spin orientation 
cotunneling can take place via both empty and doubly 
occupied dot virtual states. The corresponding domi- 
nant contribution to electric current is then proportional 
to [l/e-l/(e+U)} 2 = l/e 2 +l/(e+U) 2 -2/[e(e+U)}. It is 
thus clear that the difference in currents flowing through 
the system in the antiparallel and parallel configurations 
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FIG. 3: The spin relaxation time dependence of the lin- 
ear conductance in the antiparallel configuration (a) and the 
relative height of the zero-bias anomaly in differential con- 
ductance (b), as well as the linear-response TMR (c). The 
parameters are the same as in Fig. Q 



is equal to —2/[e(e + U)], which results from the interfer- 
ence term. Since e < and e + U > 0, this interference 
contribution is positive. As a result, the current in the 
antiparallel configuration is larger than the current in the 
parallel configuration. 

An analysis similar to that done above for differential 
conductance [see Eqs ©-©] can be performed for tunnel 
magnetoresistance. The minimum in TMR at zero bias 
in the case of a symmetric Anderson model can be then 
expressed as 



TMR,, 



2p 2 (4k B TT 2 - e 2 h/r s[ ) 
12(1 - p 2 )k B TT 2 + (3 + p 2 )e 2 h/T si ' 



(11) 
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whereas for \eV\ 3> k B T and r<l one finds 



TMR„ 



2p 2 [2(3-p 2 )k B TT 2 -e 2 h/T sf ] 

= 2(1 - p 2 )(3 - p 2 )k B TT 2 + (3 + p 2 )e 2 h/r sf ' 

(12) 

The latter formula approximates the value of TMR cor- 
responding to the bias voltage at which the differential 
conductance has a local minimum. In the slow spin re- 
laxation limit one finds TMR m j n = 2p 2 /(3 — 3p 2 ), and 
TMR max = 2p 2 /(l - p 2 ). However, in the limit of fast 
spin relaxation TMR becomes negative and is given by 

TMR min = TMR max = -2p 2 /(3 + p 2 ), (13) 

which is consistent with numerical results displayed in 

Fig.nc). 

The explicit variation of the minimum in TMR at zero 
bias with spin relaxation time is shown in Fig.nc). The 
transition between slow and fast spin relaxation in the 
dot is clearly evident. This transition is associated with 
a change of the TMR sign from positive to negative. 



B. The case of spin-split dot level 




The discussion up to now was limited to the case of de- 
generate dot level. The situation changes when e-j- 7^ ex, 
e.g., due to an external magnetic field. The magnetic field 
splits the dot level, consequently the occupation proba- 
bilities for the spin-up and spin-down electrons are not 
equal at equilibrium (V = 0). The level splitting is de- 
scribed by the parameter A = s^—s^, where the magnetic 
field is assumed to be along the magnetic moment of the 
left electrode. In Fig. 0] we show the bias voltage depen- 
dence of the differential conductance in the parallel and 
antiparallel configurations for different values of param- 
eter r. In the limit of no intrinsic spin relaxation in the 
dot (solid line in Fig. and at low bias voltage, the dot 
is occupied by a spin-up electron and the current flows 
mainly due to non-spin-flip cotunneling. The spin-flip co- 
tunneling processes are suppressed for |A| > \eV\,k B T, 
which results in the steps in differential conductance at 
A I ~ |eV|- The suppression of spin-flip inelastic cotun- 
neling was recently used as a tool to determine the spec- 
troscopic g-factor 31 . When \eV\ becomes larger than |A|, 
spin-flip cotunneling is allowed, consequently the conduc- 
tance increases. However, there is a large asymmetry of 
differential conductance in the antiparallel configuration 
with respect to the bias reversal. To understand this 
asymmetry, it is crucial to realize that when the split- 
ting A — £| — £f is larger than k B T, the single-barrier 
spin-flip cotunneling processes can occur only when the 
dot is occupied by a spin-down electron. This follows 
simply from the energy conservation rule. (The situation 
may be changed for negative A.) Thus, the single-barrier 
processes can assist the fastest double-barrier cotunnel- 
ing processes, but only for positive bias. This is because 
the fastest processes can occur when the dot is occu- 
pied by a spin-down electron for negative bias and by a 



FIG. 4: The differential conductance in the nonlinear re- 
sponse regime for different spin relaxation in the parallel 
(a) and antiparallel (b) configurations. The parameters are: 
k B T = 0.2r, e T = -16r, e x = -14T, U = 30r and p = 0.5. 



spin-up electron for positive bias, leading to larger con- 
ductance for positive than for negative bias voltage. This 
is indeed the case in the numerically calculated charac- 
teristics shown in Fig. EJb) . No such asymmetry occurs 
in the parallel configuration [see Fig. 0fa)], as now the 
system is fully symmetric with respect to bias reversal. 

The situation changes when intrinsic spin-flip relax- 
ation processes occur in the dot. For the parameters 
assumed in Fig. the spin relaxation in the dot affects 
the conductance only for \eV\ > |A|, while for \eV\ < |A| 
the conductance is basically independent of r (see Fig.0J. 
This is because for \eV\ < |A| and |A| > k B T, the dot 
is predominantly occupied by a spin-up electron and the 
transitions to the spin-down state due to relaxation pro- 
cesses are energetically forbidden. As a consequence, the 
current flows mainly due to non-spin-flip cotunneling, ir- 
respective of spin relaxation time. This scenario holds 
for both magnetic configurations of the system. 

When \eV\ > |A|, the spin-flip cotunneling processes 
can take place and the dot can be either in the spin-up or 
spin-down state. In the parallel configuration, Fig. 0Ja), 
the conductance is slightly reduced by the spin-flip relax- 
ation processes. This can be understood by realizing the 
fact that the fastest non-spin-flip cotunneling processes 
in the parallel configuration are more probable when the 
dot is occupied by a spin-down electron than by a spin- 
up one [due to smaller energy denominator, see Eq. |JI)|]. 
Since the spin-down state (as that of larger energy) re- 
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laxes relatively fast to the spin-up state (which has def- 
initely smaller energy), this leads to a reduction in the 
conductance. On the other hand, the differential conduc- 
tance in the antiparallel configuration is enhanced by the 
relaxation processes for positive bias and diminished for 
negative bias voltages. Consider first the situation for 
positive bias. As already discussed above for r = 0, an 
important role in that transport regime is played by the 
single-barrier spin- flip cotunneling processes, which open 
the system for the fast double-barrier cotunneling by re- 
versing spin of the dot from the spin-down to the spin-up 
state. The relaxation processes play a role similar to that 
of the single-barrier cotunneling, and lead to a certain in- 
crease in the conductance. For negative bias voltage, in 
turn, the fast double-barrier cotunneling processes occur 
when the dot is occupied by a spin-down electron. The 
probability of such events is decreased by spin relaxation, 
leading to a reduced conductance. An interesting conse- 
quence of the enhancement (reduction) of the differential 
conductance for positive (negative) bias voltage is an in- 
crease of the asymmetry with respect to the bias reversal 
- see the curve for r = 1 in Fig. Efb) • 



IV. QUANTUM DOTS ASYMMETRICALLY 
COUPLED TO THE LEADS 



An interesting situation occurs when the quantum dot 
is coupled asymmetrically to the left and right leads 
(Ph Pk)- For this situation we first present the re- 
sults in the case of degenerate dot level and then proceed 
with the discussion of the case when the degeneracy is 
lifted, e.g. by an external magnetic field. 



A. The case of spin-degenerate dot level 

The differential conductance for a system with one 
electrode nonmagnetic and the other one made of a fer- 
romagnet with large spin polarization (in the following 
referred to as strong ferromagnet) is shown in Fig. El 
The parts (a) and (c) correspond to the cases when 
the dot is described by an asymmetric Anderson model 
(|e| ^ e + U). The difference between (a) and (c) is due 
to different position of the dot level e, and consequently 
also e + U, with respect to the Fermi level at equilib- 
rium. More precisely, (a) and (c) are symmetric in the 
sense that |e| in (a) is equal to e + U in (c), and vice- 
versa, e in (c) is equal to e + U in (a). Part (b), in 
turn, corresponds to the dot described by a symmetric 
Anderson model, with |e| equal to e + U. 

Consider first the situation in the absence of intrinsic 
spin relaxation in the dot (solid curves in Fig.EJ. There 
is then a significant asymmetry of electric current with 
respect to the bias reversal in the situations displayed in 
Fig. Eta) and (c), whereas no such an asymmetry occurs 
for the symmetric Anderson model [part (b)]. Another 
feature of the curves shown in Fig. El is the presence of 
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FIG. 5: The differential conductance in the nonlinear re- 
sponse regime for the asymmetric (a,c) and symmetric (b) 
Anderson model for different spin relaxation. The parame- 
ters are: k B T = 0.21" , U = 60r, p L = 0.95 and p R = 0. 



zero-bias anomaly in the case described by the symmet- 
ric Anderson model [part (b)]. This anomaly is similar 
to that shown in Fig. Qta) for antiparallel configuration. 
The anomaly still exists, although it is less pronounced 
when the dot is described by the asymmetric Anderson 
model, as clearly visible in Fig. Et a ) an d (c). 

Let us discuss now the physical mechanism of the 
asymmetry with respect to the bias reversal displayed 
in Fig. Eta) and (c). When \eV\ 3> &bT, one can neglect 
the influence of single-barrier cotunnelingSS. The cotun- 
neling processes which transfer charge from one lead to 
another take place via two possible virtual states - empty 
dot (an electron residing in the dot tunnels to one of the 
leads and another electron from the second lead enters 
the dot) and doubly occupied dot (an electron of spin op- 
posite to that in the dot enters the dot and then one of the 
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two electrons leaves the dot). Consider first the situation 
shown in Fig. 0a) for positive bias (eV > 0, electrons 
flow from right to left, i.e., from normal metal to strong 
ferromagnet), and assume for clarity of discussion that 
the strong ferromagnet is a half-metallic one with full 
spin polarization (only spin-up electrons can then tunnel 
to the left lead). When a spin-down electron enters the 
dot, it has no possibility to leave the dot for a long time. 
The allowed cotunneling processes occur then via dou- 
bly occupied dot virtual states. In the absence of intrin- 
sic spin relaxation in the dot, the only processes which 
can reverse the dot spin are the single-barrier cotunneling 
ones, which however play a minor role when k^T <C |eV|. 
Thus, the current flows due to non-spin-flip cotunneling 
via doubly-occupied dot virtual states, whereas cotunnel- 
ing through empty-dot virtual states is suppressed. The 
situation is changed for negative bias (electrons flow from 
strong ferromagnet to normal metal). Now, the dot is 
mostly occupied by a spin-up electron, which suppresses 
cotunneling via doubly-occupied dot virtual state and the 
only contribution comes from cotunneling via empty-dot 
virtual state. The ratio of cotunneling rates through the 
empty dot and doubly occupied dot virtual states is ap- 
proximately equal to £ — [e/(e + U)]~ 2 . In the situation 
presented in Fig. 0a) one finds £ ^> 1. Accordingly, 
the conductance for negative bias is much larger than for 
positive bias voltage. In turn, in the case displayed in 
Fig. 0c) the ratio of the respective cotunneling rates is 
exactly equal to the inverse of this ratio corresponding to 
the case shown in Fig. 0a). Thus, now we have £ <C 1, 
and the conductance is smaller for negative bias than for 
positive. In turn, in the case presented in Fig. 0b) the 
ratio £ is equal to 1, consequently both cotunneling rates 
are equal and the conductance is symmetric with respect 
to the bias reversal. 

When \eV\ becomes of the order of k-&T or smaller, 
the rate of single-barrier cotunneling is of the order of 
the rate of double-barrier cotunneling. Therefore, the 
single-barrier processes can play an important role in 
transport. More precisely, single-barrier cotunneling pro- 
cesses can reverse spin of an electron in the dot and thus 
can open the system for the fast cotunneling processes. 
This behavior can be observed in Fig. EI a) and (c). In 
the symmetric case shown in Fig. 0b), the single-barrier 
processes lead to the zero-bias anomaly similar to that 
shown in Fig. ^a) for antiparallel configuration in a sys- 
tem with two ferromagnetic electrodes coupled symmet- 
rically to the dot. The mechanism of the anomaly is the 
same, i.e., the single-barrier cotunneling, which occurs 
for U-qT 3> |eV| , opens the system for the fast cotunnel- 
ing processes. This leads to an increase in conductance 
in the small bias range. The zero-bias anomaly exists 
also in asymmetric cases illustrated in Fig. 0a) and (c). 
However, the maximum in conductance is slightly shifted 
away from V = and is much less pronounced. 

Intrinsic spin-flip processes in the dot have similar in- 
fluence on electronic transport as in the symmetric case 
studied in the previous section. As before, relaxation 



processes remove the asymmetry with respect to the bias 
reversal and suppress the zero-bias anomaly. Thus, the 
diode-like behavior can appear only in the limit of slow 
spin relaxation, and is suppressed in the limit of fast spin 
relaxation, as shown in Fig. by the curves correspond- 
ing to r = 1. 

Some analytical expressions for the conductance can 
be obtained in the limit \eV\ 3> k^T and when assum- 
ing j»l = PR = 0. The approximate formulas for the 
differential conductance then read 



G = — — 



e 

h 2 



1 



+- 



(1- P 2 )U 
_(e + U) 2 + e 2 (e + U) 

P 2 s 2 U(e + Uf 
[e 2 {e + U) 2 - eVU 2 T 2 T si /{Ah)f 



(14) 



for positive bias, and 



G 



f!L 2 r 

h 2 
+- 



1 



(1-P 2 )U 
^\e\(e + U) 2 

p 2 \e\ z U{e + Uf 



[e 2 (e + U) 2 + eVU 2 T 2 T sf /(4h)Y 



(15) 



for negative bias. In order to describe the diode operation 
it is useful to define the ratio a of the conductances for 
positive and negative bias voltages. In the limit of long 
spin relaxation this ratio is given by 



e 2 + (1 - p 2 ){e + U)U 
(e + U) 2 + (l-p 2 )\e\U' 



(16) 



It is also worth noting that the operation of the diode 
can be tuned by the gate voltage (which shifts energy of 
the dot level). 



B. The case of spin-split dot level 

In Figure we illustrate the influence of spin splitting 
of the dot level due to an applied magnetic field, e-\ 7^ £|. 
The bias voltage dependence of the differential conduc- 
tance is shown there for different values of the parame- 
ter r, for symmetric (b) and asymmetric (a,c) Anderson 
models. The parameters used in numerical calculations 
are the same as in Fig. except for the dot level energy 
which was spin degenerate in Fig. and now is spin split. 

Let us first focus on Fig. 0a), which corresponds to 
the case of asymmetric Anderson model, e = — U /4. Two 
new features of the conductance in the limit of no intrin- 
sic spin relaxation appear in this figure. First, the con- 
ductance is now suppressed above a certain finite positive 
value of the bias voltage, and not for eV > as in the 
case shown in Fig. 0a). Strictly speaking, the conduc- 
tance drops when eV is equal to the level splitting A, 
eV « A. The shift of the main slope in the conductance 
from eV = to eV « A is due to the fact that the 
dot is blocked for cotunneling via the empty state by a 
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FIG. 6: The differential conductance in the nonlinear re- 
sponse regime for different spin relaxation. The parameters 
are: k B T = 0.2I\ e T = e - A/2, e L = e + A/2, A = 2I\ 
17 = 60r, p L = 0.95 and p R = 0. 



spin-down electron only when eV exceeds the level split- 
ting. This results from the energy conservation. Thus, 
in the case shown in Fig. EI a ) the main drop of the con- 
ductance was at eV = 0, whereas in Fig. Et a ) it is at 
eV « A. When a spin-down electron appears in the dot, 
then, for eV < A, it can always tunnel back to the source 
electrode and open the dot for fast cotunncling through 
the empty state. 

The second feature of the differential conductance in 
the r = case is the peak at eV as A, and a shallow 
plateau for \eV\ < A. The origin of the peak is similar 
to the origin of the zero-bias anomaly in the case of spin- 
degenerate dot level. On the other hand, the plateau in 
the small bias range is due to suppression of the spin-flip 
cotunncling for \eV\ < A, similarly as in Fig. How- 
ever, the plateau is now very weak and shallow. This is 
a consequence of the minor role of the spin-flip cotun- 
neling, which results from the strong spin polarization of 
ferromagnetic electrode. 



Qualitatively similar behavior may be also observed for 
e = —3U/4. This situation is shown in Fig. Et c ) and is 
analogous to the one illustrated in Fig. Etc) . In the case 
of e = —U/2 shown in Fig.Etb), the influence of magnetic 
field is qualitatively similar to that found in the case of 
a quantum dot coupled to two ferromagnetic leads in the 
antiparallel configuration, see Fig. 0] 

Transport characteristics are modified by the spin re- 
laxation in the dot. As before, the changes occur for 
\eV\ > |A|. For positive bias voltage, the spin relaxation 
processes modify the conductance mainly for e = — J7/4 
and e = -U/2 [Fig. EKb)], while for e = -311/4 
[Fig. Et c )] the changes are significantly smaller. For neg- 
ative bias, on the other hand, the conductance is slightly 
suppressed by the spin relaxation. This is because for 
negative bias the dot is mainly occupied by a spin-up elec- 
tron, and spin relaxation from the spin-up to spin-down 
dot states is suppressed due to A > IcbT. Contrarily, 
for positive bias, \eV\ > |A|, and for r s f — > oo, the dot is 
mostly occupied by a spin-down electron. Spin relaxation 
processes in the dot lead then to a certain increase in the 
occupation probability of the spin-up dot level (Pf as 1 
in the long spin relaxation time limit). This in turn may 
give rise to an increase in the differential conductance. 



V. CONCLUDING REMARKS 

We have considered numerically and analytically the 
influence of intrinsic spin relaxation in the dot on the 
differential conductance and the TMR effect in the co- 
tunncling regime. It has been shown that the zero-bias 
anomaly in the antiparallel configuration is enhanced in 
the weak spin-flip scattering regime in the dot. The op- 
timal conditions for a maximum enhancement have also 
been derived. However, in the limit of fast spin relaxation 
both the zero-bias anomaly and the dip in tunnel magne- 
toresistance at zero bias become suppressed. Moreover, 
we have found the inversion of sign of tunnel magnetore- 
sistance in the presence of fast intrinsic spin relaxation 
in the dot. 

We have also demonstrated that the zero-bias anomaly 
exists in systems with only one ferromagnetic electrode 
(the second one may be nonmagnetic). Such devices dis- 
play diode-like transport characteristics. The diode-like 
behavior exists also when the system is in an external 
magnetic field, but the operation range is changed, i.e., 
the bias voltage at which the conductance drops depends 
on the level splitting (magnetic field). Furthermore, we 
have shown that this behavior may be enhanced by in- 
trinsic spin relaxation processes in the dot. 
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